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Abstract 

The normalized differential cross section for top-quark pair production in association with at 
least one jet is studied as a function of the inverse of the invariant mass of the tf-i- 1-jet system. 
This distribution can be used for a precise determination of the top-quark mass since gluon 
radiation depends on the mass of the quarks. The experimental analysis is based on proton- 
proton collision data collected by the ATLAS detector at the LHC with a centre-of-mass 
energy of 7 TeV corresponding to an integrated luminosity of 4.6 fb”*. The selected events 
were identified using the lepton-i-jets top-quark-pair decay channel, where lepton refers to 
either an electron or a muon. The observed distribution is compared to a theoretical pre¬ 
diction at next-to-leading-order accuracy in quantum chromodynamics using the pole-mass 
scheme. With this method, the measured value of the top-quark pole mass, is: 

mP”'® - 173.7 +1.5 (stat.) +1.4 (syst.) (theory) GeV. 

This result represents the most precise measurement of the top-quark pole mass to date. 
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1 Introduction 


In the Standard Model (SM) of particle physics the couplings of the top quark to other particles are fixed 
through the gauge structure. The only free parameters in the top-quark sector of the SM are the ele¬ 
ments of the Cabibbo-Kobayashi-Maskawa mixing matrix and the top-quark mass. Due to its high value 
compared to the other quark masses, accurate knowledge of the top-quark mass is particularly relevant 
because it is related to the Higgs-boson and VT-boson masses through radiative and loop corrections. The 
precise determination of these quantities allows a stringent test of whether the model is consistent [1,2]. 
In addition, the precise knowledge of the top-quark mass is a crucial ingredient in recent evaluations of 
the stability of the electroweak vacuum [3-5]. 

The top-quark mass was determined directly at the Tevatron and at the Large Hadron Collider (LHC). A 
combination of a subset of these measurements yields a value of m, = 173.34 ± 0.76 GeV [6]. In these 
measurements, the top-quark mass is inferred from a kinematic reconstruction of the invariant mass of its 
decay products which is then calibrated to the mass definition used in the Monte Carlo (MC) simulations. 
These m, determinations lack a clear interpretation in terms of a well-defined fop-quark mass fheoref- 
ical scheme as employed in quanfum chromodynamics (QCD) perfurbafive calculafions, elecfroweak fils 
or any Iheorelical prediction in general [2, 6-9]. The values exlracled using Ihese melhods are usually 
idenlified with the top-quark pole mass but present studies estimate differences befween fhe fwo 

lop-quark mass definitions of 0(1) GeV [2,6-9]. 

The lop-quark mass can also be measured from fhe inclusive cross section for lop-quark pair {tt) pro¬ 
duction [10]. Wilh Ibis melhod fhe lop-quark mass scheme is unambiguously defined in fhe Iheorelical 
calculafions. However, lop-quark mass delerminafions based on cross-seclion measuremenfs are less pre¬ 
cise, in fheir currenf form, fhan the other techniques based on kinematic reconstruction. This is due to a 
relatively weak sensitivity of the inclusive top-quark pair production cross section to the top-quark mass, 
as well as to the large uncertainties on the factorization and renormalization scales and the proton parton 
distribution function (PDF). To date, the most precise measurement of this type is based on the 7 TeV and 
8 TeV data samples collected by the ATLAS experiment during the years 2011 and 2012, which yields 
^poie _ j' 72.9]]26 GeV [11]. The results from the CMS experiment using this technique only include data 
collected during 2011 [12,13]. 

In this paper the method described in ref. [14] is followed. The top-quark mass is extracted from a 
measurement of the normalized differential cross section 'R(m[’°*^,pi) for tt production with at least one 
additional jet, ti + 1-jet, as a function of the inverse of the invariant mass of the ti + 1-jet system, ps <x 
1 / Vvr+i-jet- This distribution is sensitive to the top-quark mass because the amount of gluon radiation 
depends on its value, with large effects in the phase-space region relatively close to the tf-i- 1-jet production 
threshold. This method combines the rigorous interpretation of the mass inferred from the inclusive cross 
section with the advantage of a greater sensitivity. 

The measurement is performed using 7 TeV proton-proton (pp) collision data collected by the AT¬ 
LAS experiment [15], corresponding to an integrated luminosity of 4.6 fb“^ with a total uncertainty of 
±1.8% [16]. Events are selected by using the “lepton±jets” final sfafe fo identify the tt system and at least 
one additional jet. In this channel one W from the top decay produces a lepton (electron or muon) and 
a neutrino whereas the other W produces a pair of light quarks. Events are thus required to have exactly 
one lepton, two jets identified as f^-quark jefs, af leasf fhree addifional jefs nof idenfified as f7-quark jefs 
and a significanf amounf of missing fransverse momenfum due fo fhe neufrino fhaf escapes defec- 

fion. The differenfial cross section is correcfed fo fhe parfon level (a procedure referred fo as unfolding 
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in the following) and normalized. The top-quark pole mass is then extracted through a fit to the data us¬ 
ing the predicted normalized differential tt + 1-jet cross section from a next-to-leading-order calculation 
combined with parton showering (NLO-tPS) [14, 17-19]. 


2 Definition of the observable 

The method to extract the top-quark pole mass followed here, proposed in ref. [14], exploits the fact that 
the top-quark mass dependence of the tt + 1-jet cross section, crtf+i-jeu is enhanced in the phase-space 
region relatively close to the ti+ 1-jet production threshold. This method uses the predictions for tt+ 1-jet 
production at hadron colliders at NLO accuracy reported in refs. [17,18]. A well-defined top-quark pole 
mass can be extracted by comparing these calculations with the measurement of the normalized ti + 1-jet 
cross section in pp collisions as a function of the inverse of the invariant mass of tho ti + 1-jot 

system: 


where ps is defined as 




0"«r+i-jet dp 


dcTf^+i.jet pole , 


Ps = 


2mo 


( 1 ) 

( 2 ) 


and niQ is an arbitrary constant of the order of the top-quark mass. Here and in the following, mq = 
170 GeV is used. The anti-kj jet reconstruction algorithm [20,21] is employed to reconstruct the jets. The 
extra jet, beyond those which originated from the ti decay, is the leading jet with a transverse momentum 
Pt > 50 GeV and a pseudorapidity \p\ < 2.5. 1 The observable K defined in this way is infra-red safe as 
demonstrated in the study of ref. [14]. 

In this analysis, the measured normalized and differential cross section, unfolded to parton level, is 
compared to the theoretical calculations at NLO accuracy, after adding the parton shower evolution 
(NLO-i-PS). Including the parton shower is expected to give a better description of the final-state phase 
space than the NLO calculation alone and is implemented in the publicly available MC generator de¬ 
veloped in ref. [19]. This generator uses Powheg (PowHEG-ttJ) [19,22,23] matched with the Pythia 
v8 [24] parton shower. Using a fixed order NLO calculation to fit the data gives a similar ^^-distribution 
but leads to an estimated top quark pole mass about 0.3 GeV lower than using a NLO-i-PS calculation. 

This difference is well below the present theoretical uncertainty of the calculation. Differences due to the 
use of Pythia v8 or Pythia v6 [25] are below this value of 0.3 GeV. 

In the NLO calculation, it is assumed that the top quarks are stable. Possible effects due to radiation from 
top-quark decay products and virtual corrections to the decay are small compared to the overall theoretical 
uncertainty. Quantum chromodynamics corrections to the decay do not affect the mass renormalization of 
the top quark at the same order of accuracy as considered in the calculation because the renormalization is 


* ATLAS uses a right-handed coordinate system with its origin at the nominal interaction point (IP) in the centre of the 
detector and the z-axis along the beam pipe. The x-axis points from the IP to the centre of the LHC ring, and the ^-axis points 
upward. Cylindrical coordinates (r, are used in the transverse plane, 0 being the azimuthal angle around the beam pipe. The 
pseudorapidity is defined in terms of the polar angle 6 as ?/ = - In [tan (6/2)]. Transverse momentum and energy are defined 
as pt = P sin 6 and Ej = E sin 6, respectively. 
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purely determined from the QCD self-energy corrections of the top-quark propagator, which is included 
in the calculation. Furthermore, recent calculations in refs. [26, 27] include NLO QCD corrections to 
the total and differential tt cross section assuming the top quarks to be off-shell. In this framework, the 
results in the on-shell approximation are reliable and off-shell effects are at the sub-percent level. In 
the following, it is assumed that similar results hold for the quantity H. In fact, as H is expressed as a 
normalized cross section, even smaller effects can be expected due to possible cancellations. 


3 The ATLAS experiment 

The ATLAS experiment [15] is a multipurpose detector with a forward-backward symmetric cylindrical 
geometry. The inner tracking detector (ID) consists of a silicon pixel detector, a silicon microstrip detector 
(SCT), and a straw-tube transition radiation tracker (TRT). The ID is surrounded by a thin superconduct¬ 
ing solenoid which provides a 2 T magnetic field and by a high-granularity liquid-argon (LAr) sampling 
electromagnetic calorimeter. The electromagnetic calorimeter is divided into a central barrel (|? 7 | < 1.475) 
and end-cap regions at each end of the barrel (1.375 < \t]\ < 2.5 for the outer wheel and 2.5 < \r]\ < 3.2 for 
the inner wheel). A steel/scintillator-tile calorimeter completes the measurement of hadronic showers in 
the central pseudorapidity range (It/I < 1.7), while a LAr hadronic end-cap calorimeter provides coverage 
over 1.5 < \t]\ < 3.2. The forward regions (3.2 < \t]\ < 4.9) are instrumented with LAr calorimeters for 
electromagnetic and hadronic measurements. The muon spectrometer (MS) surrounds the calorimeters 
and consists of three large air-core superconducting toroid systems providing bending powers of 3 Tm 
in the barrel and 6 Tm in the end-caps, a system of precision tracking chambers (|?/| < 2.7), and fast 
detectors for triggering (\t]\ < 2.4). The combination of all these sub-detectors provides charged-particle 
measurements together with efficient and precise identification of leptons and photons in the pseudorapid¬ 
ity range of \t]\ < 2.5. Jets and L"™®* are reconstructed using energy deposits over the full coverage of the 
calorimeters, \r]\ < 4.9. The reconstructed muon momenta are also included in the evaluation of 
Evaluation of the luminosity scale is performed using several luminosity-sensitive detectors. The ATLAS 
experiment has a three-level trigger system [15]. The first-level trigger is hardware-based and uses a sub¬ 
set of the detector information to reduce the event rate to at most 75 kHz. The second and third levels are 
software-based and together reduce the event rate to about 300 Hz. 


4 Data sample and Monte Carlo simulation 

The data considered in this analysis correspond to an integrated luminosity of 4.6 fb“^ of proton-proton 
collisions at a centre-of-mass energy of 7 TeV. They were recorded in 2011 during periods with stable 
beam conditions and with all relevant subdetector systems operational. The events were selected by 
single-lepton triggers that require a minimum transverse momentum of 18 GeV for muons and a minimum 
of 20 to 22 GeV for electrons, depending on the data-taking conditions. 

In this analysis, several MC samples are used for the modelling of tt pair production and the main back¬ 
ground processes. For the simulation of the tt signal and the unfolding, the Powheg code (PowuEG-hvq, 
patch4 [28]) is used to calculate the QCD matrix element at NLO with the CTIO [29] PDF set. The 
parton shower and the underlying event are added using the Pythia v6.4 [25] generator with the Perugia 
2011C parameter set (tune) [30]. Several MC sets of events were generated using different top-quark 
masses. The nominal MC sample, which is used in the present study to compare the MC predictions 
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with data, is the largest and is produced assuming a top-quark mass of nit = 1V2.5 GeV. The corres¬ 
ponding cross section of the nominal tt sample is 177^|j (theo.) pb as predicted by the calculations in 
refs. [31-35], which include the next-to-next-to-leading order (NNLO) and the resummation of next-to- 
next-to-leading logarithmic (NNLL) soft gluon terms with top-i--i-2.0 [34]. In addition to this nominal 
sample, five other samples are employed with different assumptions about the input top-quark mass in the 
range from 167.5 GeV to 180 GeV in steps of 2.5 GeV. 

Background processes to the tt + 1-jet final state under study are simulated with various MC generat¬ 
ors. Single top quark production in the s-, Wt- and t-channels is simulated using Powheg matched with 
Pythia v6.4. The Perugia 2011C tune is used. The production of W/Z bosons in association with jets 
(VP-tjets or Z-i-jets) is simulated using the ALPGEN generator (v2.13) [36] with the leading-order (LO) 
CTEQ6L1 [37] PDE set. These calculations are interfaced with Herwig 6 [38] for the parton shower 
and Jimmy v4.31 [39] for the underlying-event modelling. IT-i-jets events containing heavy-flavour quarks 
are generated separately using leading-order matrix elements with massive b- and c-quarks. Possible 
double-counting due to heavy quarks produced by the parton shower is considered and removed. The 
total number of IT-i-jets events is normalized by exploiting the lepton charge asymmetry observed in data, 
following the method in ref. [40]. Diboson events {WW, ZZ, WZ) are generated using Herwig 6 with the 
MRSTMCal [41] PDE The background from misidentified and non-prompt leptons is estimated using a 
data-driven matrix method described in ref. [42]. 

Multiple soft pp interactions generated with Pythia v.6.425 using the AMBT2B tune [43] are added to all 
simulated events in order to account for the effect of multiple pp interactions in the same and neighbouring 
bunch crossings (pile-up). 

The response of the ATEAS detector is simulated using a detailed description of the detector geo¬ 
metry [44] in GEANT4 [45]. Simulated events are reconstructed using the same software as used for 
the data. 


5 Object definition and basic selection 

The analysis applies several requirements to the events and makes use of reconstructed electrons, muons, 
jets and E™®®. Electron candidates are reconstructed from energy deposits in the electromagnetic calori¬ 
meter using criteria based on the shower shape, and they must be matched to a charged-particle track in 
the ID [46]. Electrons must have a transverse momentum of px > 25 GeV and I 77 I < 2.47. Events with 
electrons falling in the calorimeter barrel/end-cap transition region, corresponding to 1.37 < |p| < 1.52, 
are rejected. Muon candidates are identified by combining track candidates in the MS and the ID [47]. 
All muons are required to have a transverse momentum pj > 25 GeV and \r]\ < 2.5. All muons must 
additionally satisfy a series of selection criteria on the number of hits per track in the various tracking 
sub-detectors [48]. Einally, electrons and muons have to match corresponding objects that have fired the 
trigger. 

Isolation criteria are applied to electron and muon candidates to reduce the background from hadrons 
mimicking lepton signatures and backgrounds from heavy-flavour decays inside jets [49]. 

Energy deposits in the calorimeters are combined into three-dimensional clusters [50]. Erom these 
clusters, jets are reconstructed using the anti-k^ jet algorithm with a radius parameter of 0.4. Recon¬ 
structed jet energies in simulations are calibrated from stable-particle jets. Residual calibrations, derived 
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by using in situ methods where the jet transverse momentum is compared to that of a reference object 
(e.g. using y/Z+jet events), are then applied to data relative to the simulations [51]. 

Reconstructed jets must have px > 25 GeV and \r]\ < 2.5. To suppress the contribution from low-px 
jets originating from pile-up interactions, a jet vertex fraction requirement is applied: at least 75% of the 
summed scalar pj of tracks associated with the jet must be due to tracks originating at the event primary 
vertex. This primary vertex is defined as the vertex with the highest ZtrkCPx^)^ among all candidates with 
at least five associafed fracks (frk) wifh p'^ > 0.4 GeV [52]. Jefs confaining fj-hadrons are idenfified 
using a neural nefwork exploifing fhe long lifelime of fj-hadrons af a working poinf resulfing in a fagging 
efficiency of 70% in simulafed tt evenfs [53-55]. 

The fransverse momenfum of neufrinos escaping fhe defector is assumed to be identical fo which is 
reconsfrucfed as fhe magnifude of fhe momenfum imbalance in fhe fransverse plane as obfained from fhe 
negafive vector sum of fhe momenfa of all energy deposifs. If is reconsfrucfed from fopological clusfers, 
calibrafed for elecfromagnefic objecfs and correcfed according fo fhe energy scale of fhe idenfified objecfs. 
Muons confribufions are also included by using fheir momenfum measured in fhe inner defecfor and fhe 
muon specfomefer [56]. 

Evenfs musf nof confain jefs wifh px greafer fhan 20 GeV arising from ouf-of-fime energy deposifs or 
from energy deposifs wifh a hardware or calibration problem. Exacfly one isolafed elecfron or muon, and 
af leasf five jefs are required wifh exacfly fwo of fhe jefs fagged as f?-quark jefs. The magnifude of fhe 
missing fransverse momenfum and fhe fransverse mass of fhe sysfem formed by fhe charged lepfon and 
fhe neufrino,^ ni^, musf bofh exceed 30 GeV. 


6 Reconstruction of the + 1-jet system 


Affer fhe basic selection, a kinemafic reconsfrucfion of fhe evenfs is performed fo idenfify fhe IT-boson 
and fop-quark candidafes. The lepfonically-decaying W boson is idenfified wifh fhe charged lepfon and 
fhe neufrino, where fhe longifudinal momenfum is inferred using a consfrainf on fhe IT-boson mass. 
Candidafes for fhe hadronically-decaying W boson are consfrucfed by considering all possible pairs of 
jefs among fhose nof idenfified as f^-jefs. Accepfed evenfs musf fulfil fhe following condifions: 


m 


ref 

W 


0.9 < a = — < 1.25, 


ni; 


(3) 


Akt{i,j) = min{pj, p!^) ■ AR{i,j) < 90 GeV, (4) 

where fhe indices i and j run over all jefs nof idenfified as f^-jefs, is fhe invarianf mass of jefs i and 
j, and = 80.4 GeV [57]. All combinations of fhe fwo f^-jefs wifh fhe VT-boson candidafes are con¬ 
sidered as fop-quark candidafes. Once bofh requiremenfs are applied, fhe permufafion fhaf minimizes fhe 
invarianf mass difference of fhe hadronically and lepfonically decaying fop-quark candidafes is selecfed. 
The applicafion of fhese fwo requiremenfs reduces fhe combinatorial background. 


^ The VT-boson transverse mass is defined as px.vll - cos — (py)] where t and v refer to the selected lepton and 

the neutrino whose information is associated with the vector, respectively. 


6 





Events Uncertainty 

Signal (tt, nit - 172.5 GeV) 

2050 ± 320 

VT-i-jets 

Z-i-jets 

Single top (m^ = 172.5 GeV) 

ITW, ZZ, IVZ 

Misidentified and non-prompt leptons 

31 ± 16 

6 ± 4 

62 ± 34 

1 ± 1 

22 ± 13 

Total Background 

Total Predicted 

Data 

121 ± 40 

2170 ± 320 

2256 


Table 1: Event yields and their uncertainties after the reconstruction of the tt+ 1-jet system. The quoted uncertainty 
values include the total statistical and systematic uncertainties as described in section 8. 







150 200 250 

b-tagged jets Pj [GeV] 




Figure 1; The data for various kinematic distributions (transverse momentum, pj, of the lepton, pj of all the b- 
tagged jets, pj of all non-Zj-tagged jets and the total jet multiplicity) are compared to the nominal tt MC sample 
(Powheg-i-Pythia) plus backgrounds after the final kinematic reconstruction of the tl + 1-jet events. The total 
background estimated in table 1 is shown in dark grey. The uncertainty band includes the total statistical and 
systematic uncertainties as described in section 8. 
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Figure 2; The data for various kinematic distributions (the reconstructed mass of the hadronically and leptonically 
decaying top-quark candidates, the pj of the additional jet and the invariant mass of the tt -n 1-jet system) are 
compared to the nominal ff MC sample (Powhegh-Pythia) plus backgrounds after the final kinematic reconstruction 
of the tt -I- 1-jet events. The total background estimated in table 1 is shown in dark grey. The uncertainty band 
includes the total statistical and systematic uncertainties as described in section 8. 


The magnitude of the momentum vectors of the light-quark jets associated with the hadronic W boson 
are scaled using the value of a derived from eq. 3. A further requirement on the ratio of the reconstructed 
top-quark kinematic mass for the leptonic and hadronic decays, > 0.9, is imposed to 

increase the signal to background ratio. The leading-pj jet is selected among the remaining jets and is 
identified as the extra jet completing the tt + 1-jet system. The extra jet must satisfy pj > 50 GeV and 
|77| < 2.5. 

The event yields after the final selecfion are presented in fable 1 . The number of selecfed dafa evenfs is in 
good agreemenf wifh fhe MC expecfafion. 

Afler applying all fhe selecfion criteria, data were compared to the expectations. A representative subset 
of kinematic and angular distributions is shown in figures 1 and 2. In all cases, good agreemenf befween 
dafa and predicfion is observed. In addition fo fhese plofs, fhe number of reconsfrucfed tt+ I-jef evenfs as 
a function of is shown in figure 3. 






















Figure 3: Number of reconstructed events as a function of (niQ - 170 GeV) related to the inverse of the invariant 
mass of the tt + 1-jet system. The data are compared to the nominal ff MC prediction (Powheg-i-Pythia) plus 
backgrounds, which assumes a top-quark mass m, =172.5 GeV after the final kinematic reconstruction of the 
tt + 1-jet events. The backgrounds and the systematic uncertainties are estimated as in figure 1 and 2. 


7 Top-quark mass determination 

The top-quark pole mass is extracted by fitting the parameterized NLO-i-PS prediction to the measured 
distribution of the normalized differential cross section K defined in eq. 1 affer background subfracfion 
and correcfion for defector effecfs and hadronizafion. This mefhod follows a similar procedure to fhaf 
employed in ref. [58]. The observed number of tt + l-jef evenfs in figure 3 is used to consfrucf fhis 
‘R-disfribufion af parfon level. 

The recorded parfon-level information in fhe nominal ATLAS signal MC sample (see section 4) did 
nof allow fhe consfrucfion of fhe exfra jef af fhe parfon level, as required by fhe fheorefical calculafion. 
Consequenfly, a simple direcf connection befween fhe ti + l-jef sysfem af defecfor level and af parfon 
level could nof be made. For fhaf reason, an intermediate sfafe corresponding fo fhe firsf gluon emission 
af parton level {tt -i-g) was infroduced to bridge fhe connection and an additional MC sample of evenfs 
was generated. This laffer sample Iransformed fhe 'R-disfribufion of fhe tt+g sysfem info fhe final K- 
disfribufion of fhe tt + l-jef sysfem, which fhis time was defined as in fhe fheorefical calculafion. This 
sample was validafed using ATLAS procedures and did nof include fhe full defecfor simulafion. The final 
correction procedure confained fhe following sfeps: 

tt + l-jef (defecfor level) ^ tt+g (parfon level) tt + l-jef (parfon level). 

As a firsf sfep fhe migrafion mafrix relafing fhe disfribufion of tt + l-jef evenfs af defecfor level and fhe 
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tt +g events at parton level was calculated using the nominal MC simulation. Data were grouped in 
bins as a function of ps with a variable bin size as displayed in table 2 and in figure 4. This choice was 
the result of a compromise between having values of the diagonal terms in the migration matrix above 
50% and optimizing the sensitivity of the “^^-distribution to the top-quark mass according to the study of 
ref. [14]. The migration matrix was then inverted using the Singular Value Decomposition (SVD) method 
described in ref. [59]. This algorithm minimizes the statistical fluctuations inherent in the matrix-inversion 
process. 

The above correction restricts the kinematical phase-space region to that of the events satisfying the 
selection criteria. Hence an additional correction is needed to extend this region to the full acceptance 
considered by the theoretical calculation (acceptance term). Finally an additional step was implemented 
to convert the “R-distribution corresponding to the tt+g system into the “R-distribution of the tt + 1-jet 
system defined as in the theoretical calculation. The final unfolding procedure is described as follows: 

)) 

where ®is the term describing the transformation of the tt + 1-jet distributions from 
detector level to the parton level at its first gluon emission {tt+g). The acceptance term is 

J ’ 

and deviates from unity by less than 2% over the full ps range. The factor 

transforms the ^^-distribution which corresponds to the tt+g system into the ^^-distribution of the tt+ 1-jet 
system. This correction factor is typically within 10% of unity. 

In more detail, the first part of eq. 5 corresponds to: 


^cor-data(^^) = 




,^det-data^^ 


Clips) 

n’’^^iPs) 


n^^^iPs) 


(5) 


(Af-^ 


1 


ZM-. 

J 


1 , 


J 


N 


ff+l-jet 


^P’s 


( 8 ) 


where i,j refers to the bin numeration defined for the ps variable, is the number of tt + 1-jet 

events reconstructed (background subtracted) in the y-th bin, number of reconstructed 

tt + 1-jet events (background subtracted), Mt) is the inverse of the migration matrix and Ap^ is the width 
of the /-th bin. 


Once data were properly corrected, the value of was determined by fitting the unfolded ^^-distribution 

with the NLO-i-PS prediction from ref. [19] using the least-squares method. Table 2 shows the predicted 
values of *R (y?*eory^ diffej-gut values. The fit minimized ax^ defined as: 


= 2^(y?cor-data _ y^theory^^pole^^^-l^^cor-data _ ^theory 

u 


(9) 
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(^theory 

pole 

Ps {mo = 17()73eV7~~-------,^ 

170 GeV 

172.5 GeV 

175 GeV 

177.5 GeV 

180 GeV 

0 to 0.25 

0.1327(9) 

0.1390(44) 

0.1425(8) 

0.1487(10) 

0.1548(6) 

0.25 to 0.325 

1.104(14) 

1.134(6) 

1.172(13) 

1.208(15) 

1.251(9) 

0.325 to 0.425 

1.972(9) 

2.027(4) 

2.070(9) 

2.130(11) 

2.185(7) 

0.425 to 0.525 

2.506(12) 

2.561(6) 

2.587(11) 

2.644(13) 

2.674(8) 

0.525 to 0.675 

2.143(8) 

2.125(4) 

2.116(7) 

2.085(9) 

2.060(6) 

0.675 to 1.0 

0.353(2) 

0.316(1) 

0.287(2) 

0.252(2) 

0.223(1) 


Table 2: The ‘??-distribution calculated using generated ti + 1-jet samples at NLO-i-PS accuracy for different 
values at parton level (corresponding to gq gy ^j^g q^oted uncertainties in parentheses reflect the 

statistical precision of the calculation. 


where is the data value in the /-th bin of the corrected ‘/^-distribution and V ^ is the inverse of 

the statistical covariance matrix of the unfolded '^-distribution. The quantity represents the 

theoretical prediction for the /-th bin and contains the dependence on the top-quark pole mass. The covari¬ 
ance matrix is obtained by producing a sample of 500 pseudo-experiments scattered around the measured 
values of the ^^-distributions assuming Gaussian statistical errors. Each of these pseudo-'R-distributions 
was then corrected following the unfolding procedure described above and finally the covariance matrix 
(V) was evaluated accordingly. 

The inferred value is the one which minimizes the in eq- (9) calculated by considering all bins 
except for the least sensitive one (0 < < 0.25). This is done because the ^^-distribution is constrained 

by the normalization condition. The extracted mass value does not significantly depend on this bin choice. 
The selected configuration is the one which gives the highest expected precision. The statistical uncer¬ 
tainty is taken as the mass shift that increases the;y^ by one unit with respect to the minimum (Ax^ = +1)- 
The possible impact of non-perturbative effects near the threshold for tt + 1-jet production (p^ ~ 1) was 
also studied by choosing a restricted range of ps from 0 to 0.9 and by studying the results for each bin 
independently. No significant effect is observed. 

The corrected ‘/^-distribution is shown in figure 4. For comparison purposes the predictions for three 
top-quark pole mass values are also shown -170, 175 and 180 GeV) together with the top-quark 
mass extracted from the best fit, which is =173.7±1.5 (stat.) GeV. Only the statistical uncertainty is 
shown in this figure. 

The corrections for detector and hadronization effects based on the MC simulation might introduce a 
dependence on the input top-quark mass assumed in the generator. The possible impact was quantified 
by generating input tt+ 1-jet distributions using fully simulated samples with Powheg -i- Pythia and with 
different MC masses ranging from 167.5 GeV to 180 GeV, i.e. keeping the same nominal MC parameter 
set except for the input top-quark mass. Each of the corresponding distributions is unfolded with the 
same procedure as used for data, fixing the migration matrix and correction factors, which are defined for 
a fixed top-quark mass of m, = 172.5 GeV. Fits to the resulting ^^-distributions are performed using the 
parameterised mass dependence of the theoretical predictions obtained from the MC samples with input 
top-quark masses between 167.5 GeV and 180 GeV. Each top-quark mass extracted is compared with the 
corresponding input top-quark mass. Figure 5 shows the difference between the input and fitted masses as 
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Figure 4; 'R-distribution at parton level corrected for detector and hadronization effects after the background sub¬ 
traction as a function of ps (/«o = 170 GeV). The predictions of the tt + 1-jet calculation at NLOh-PS using three 
different masses =170, 175 and 180 GeV) are shown together with the result of the best fit to the data, 

= 173.7+1.5 (stat.) GeV. The black points correspond to the data. In the lower part of the figure, the ratios of the 
different '^-distributions to the one corresponding to the best ht are shown. The shaded area indicates the statistical 
uncertainty. 


a function of the input mass. In the range studied here, all fit results are compatible with the input values 
within their statistical uncertainties. 

Existing generated samples with Powheg + Herwig 6 including full ATLAS simulation were used to 
make the correction of the data without using the intermediate state of the tt+g system. This cross¬ 
check allowed investigations of potential biases introduced by this step. When using this sample the 
correction procedure was tested including and excluding the tt+g intermediate state. The two methods 
gave compatible results within ~0.1 GeV for well within the statistical precision of the test, ~0.25 
GeV. 
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Figure 5; Difference between the fitted mass and the top-quark mass assumed in the generated ff MC predictions 
(Powheg-hPythia) including full detector simulation as a function of the input mass. The same unfolding procedure 
employed for data is performed in this study. The migration matrix and correction factors are defined for a fixed 
top-quark mass of m, - 172.5 GeV. The fit is performed using the parameterised mass dependence of the theoretical 
predictions obtained from the MC samples. A fit to a straight line including the point at 172.5 GeV is performed. 
The obtained mean value and;t'^/NDF are shown. 


8 Statistical and systematic uncertainties 


This section describes the uncertainties that affect the extraction of the top-quark pole mass. The stat¬ 
istical uncertainty of the corrected result is evaluated using eq. 9 using toy MC experiments to derive 
the covariance matrix of the fit. The additional (small) uncertainty due to the limited number of MC 
events used to define the unfolding procedure is evaluated by varying the migration matrix according to 
their statistical uncertainties. The systematic uncertainties are split into four categories: theoretical un¬ 
certainties, signal- and detector-modelling uncertainties, and finally background uncertainties. They are 
described in the following subsections and summarized in table 3. 


8.1 Theoretical uncertainties 

Scale variations: The calculation of K is performed by setting the renormalization scale Qxr) equal to 
the factorization scale (/ip)- To estimate the uncertainty due to the missing higher-order terms in the 
calculation, these scales are varied around the central values ft - by a factor of two 

up and down. The data were fitted using the predictions with a scale twice or half the nominal 
value (yu = // = The alternative choices for the scale lead to a 0.44 GeV lower 

value for the top-quark pole mass for ft - and to a 0.93 GeV higher value for ft - jl. 
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In principle one could also vary the factorization scale independently from the renormalization 
scale. This exercise was considered in ref. [19] and the results obtained showed very good agree¬ 
ment with those from the restricted scale variation considered in ref. [18]. As a consequence no 
significant changes in the estimation of this uncertainty are expected when considering an inde¬ 
pendent variation of the scales. 

It is often argued that for normalized cross sections the method described above to evaluate the 
effect of uncalculated higher order terms in the perturbative calculations might be unrealistic and 
reduce its dependence due to cancellations of ag in the ratio. Therefore a different approach to 
evaluate this uncertainty was considered in ref. [14]. It consisted of expanding the '^-distribution in 
powers of as, thus avoiding the ratio. It was found that the two methods gave consistent estimates 
of the uncertainty. 

In addition to these cross-checks the size of the NLO correction with respect to the LO was also 
computed for "R. The study compared the LO-tPS prediction using a fixed scale jj. = mt and a 
variable dynamic scale // = ^ fo fhe NLO-i-PS predicfion wifh fixed scale. Differences 

in fhe range from 0.6 fo 0.8 GeV were observed. The small NLO correction indicafes fhaf fhe 
calculafion converges well. 

Proton PDF and agi The uncerfainfies on fhe proton PDF and on fhe value of fhe sfrong coupling con- 
sfanf as used in fhe tl + 1-jef calculafion are propagafed by tiffing fo various 'R-disfribufions im- 
plemenfed in NLO-i-PS calculations using differenf PDF sefs wifh differenf as values. The cenfral 
CTIO [60], MSTW2008nlo90cl [61,62] and NNPDF [63] PDF sefs were employed fo esfimafe 
this uncertainty. For each of these sets the central value of the resulting top-quark mass was cal¬ 
culated and the uncertainty due to the PDF corresponds to half of the maximum difference. The 
impacf of varying only as was esfimafed and found fo be very small, = 0.01 GeV (for 

Aas = ±0.002 and fhe CTIO PDF sef) since fhe dependence of R on as nearly vanishes due fo fhe 
use of a normalized differenfial disfribufion. 

The fofal fheorefical uncerfainfy is fhe sum in quadrafure of fhe confribufions of scale variafions and fhe 

PDF and as uncerfainfies. 


8.2 Detector modelling 

The uncertainties on the reconstruction efficiency and the energy measurement of basic reconstructed 
objects (leptons, and jets) are propagated to the uncertainty on the value of the top-quark mass. 
Variations of all these quantities by ±1 standard deviation are implemented in MC samples that are then 
unfolded using the nominal response matrix. A fit to the resulting ^^-distribution is performed and the 
top-quark mass is extracted. In the following — unless otherwise stated — the systematic uncertainties 
arising from the different modelling sources are calculated as half of the difference between the upward 
and downward variations. 

Jet energy scale (JES) and A-jet energy scale: To estimate the impact of the jet energy scale uncer¬ 
tainty on the result, the jet energy is scaled up and down by its uncertainty for 21 uncorrelated 
components which are considered separately [51,64]. 

These are the experimental sources of uncertainty with the largest impact on the precision of the 
mass determination. 
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Jet energy resolution and jet reconstruction efficiency: The effect of the jet energy resolution uncer¬ 
tainty is evaluated by smearing, before the event selection, the energy of the jets by a Gaussian 
function with a width chosen in agreement with the jet energy resolution uncertainty. The effect of 
the jet reconstruction ethciency uncertainty is evaluated by randomly discarding a fraction of jets 
from the events before the selection [65]. 

ft-tagging efficiency and mistag rate: Differences in the Zr-tagging efficiency as well as c-jet and light- 
jet mistag rates in data and simulation are parameterized using correction factors, which are func¬ 
tions of pt and rj. These corrections are derived from data including tt events and they are varied 
by their uncertainties (see refs. [53-55]). Similarly to the JES uncertainty, the uncertainty on the 
correction factors for the t7-tagging efficiency is separated into several uncorrelated components. 
The systematic uncertainty is assessed by changing the correction factor central values by ±lo- for 
each component, and performing the mass extraction. The final uncertainty due to the Z^-tagging 
efficiency is calculated as the quadratic sum of all contributions. 

Lepton identification and lepton energy resolution: The correction factors applied to the lepton iden¬ 
tification are measured by comparing high-purity events from simulation and data including Z, W 
and //i/r decays for electrons [46], and Z, W, Jlijf and Y decays for muons [48]. For the meas¬ 
urement of the lepton energy or momentum scale uncertainties, a similar procedure is used. The 
systematic uncertainties associated to each correction factor are estimated by changing its central 
value by ±lo- and applying the mass extraction procedure. The quadratic sum of all contributions 
gives the quoted final uncertainty due to this source. 

Modelling of the Uncertainties on the energy scale of jets or leptons are also propagated to the 

uncertainty of the Other contributions to this uncertainty originate from the energy scale and 
resolution of the soft calorimeter energy deposits which are not included in the reconstructed jets 
and leptons, and contribute only to the estimation of The is scaled up and down by its 
uncertainty. The effect of these changes is propagated to the simulated events allowing to evaluate 
the impact on the top-quark mass measurement following the same procedure as for the rest of 
systematic uncertainties. 


8.3 Signal modelling 

The signal modelling uncertainties originate from: the choice of matrix element, the parton shower and 
hadronization model, and the choice of the PDF set used in the simulation of tt events. In addition, uncer¬ 
tainties on the modelling of the initial- and final-state QCD radiation (ISR/FSR), of colour reconnection, 
and of the underlying event are also accounted for. Their impact on the extracted mass is estimated using 
alternative MC samples. The alternative ^^-distribution samples are corrected using the nominal response 
matrix and the deviation from the result of the nominal MC sample is used to estimate the uncertainty. 

MC generator and hadronization: The uncertainty associated with the choice of MC generator is eval¬ 
uated by comparing two NLO MC generators interfaced to the same parton shower and hadroniza¬ 
tion program: Powheg-box [22,23] and MC@NFO [66] both interfaced to Herwig 6 are compared. 
The difference between the extracted masses is taken as the generator uncertainty. 

The uncertainty associated with the hadronization is estimated by comparing the results obtained 
with PowHEG interfaced to either Pythia or Herwig 6. The full difference is quoted as the hadron¬ 
ization uncertainty. 
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Initial- and final-state radiation (ISR/FSR): The effect of the ISR and FSR modelling uncertainties is 
evaluated by comparing two simulated signal samples with varied radiation settings. The samples 
to evaluate the ISR/FSR uncertainty are generated with Alpgen(v2.13) [36]-i-Pythia, which is a 
multileg MC generator that generates, at LO, tt plus up to five partons. The samples to estimate 
the ISR/FSR uncertainty correspond to variations of the KTFAC parameter in Alpgen between a 
factor two up and down of its nominal value (with the Perugia 2011 radLow and radHi tunes 
respectively [30]). This parameter determines the scale at which as is evaluated for additional 
gluon emissions and the size of the variation considered is compatible with the measurements of 
additional jet activity in tt events [67]. The ISR/FSR uncertainty is evaluated by taking half the 
difference between the fitted top-quark masses from the two samples. 

Colour reconnection and underlying event: The impact of the uncertainties in the MC models describ¬ 
ing colour reconnection and the underlying event is estimated by comparing several Powheg MC 
samples with different tunes. The effect of the colour reconnection modelling uncertainty is estim¬ 
ated as the difference between the result obtained with the nominal Powheg sample with the Perugia 
2012 (P2012) tune and an alternative sample with the Perugia 2012 loCR tune [30]. To estimate 
the uncertainty on the underlying event modelling, the Perugia 2012 mpiHi tune [68] is compared 
with the P2012 tune. In both cases the full mass difference from the default value is taken as the 
systematic uncertainty. 

Proton PDF : Uncertainties on the proton PDF give rise to uncertainties on the efficiency of the basic 
event selection. These uncertainties are calculated following the PDF4LHC recommendations [69] 
using tt events simulated by MC@NLO interfaced to Herwig 6. This uncertainty accounts for 
the effects of the PDF on the theoretical modelling of the tt system, the hadronization and the 
experimental data analysis. To a large extent the first contribution is already considered in the 
evaluation of the theoretical uncertainties of section 8.1. In the present work the two uncertainties 
are considered independently. This evaluation of the total PDF uncertainty is therefore regarded as 
a conservative approach. Using the values from table 3 and considering both uncertainties either 
completely correlated or uncorrelated changes the overall uncertainty on the PDF from 0.54 GeV 
to 0.58 GeV which has a rather minimal impact. 


8.4 Background modelling 

The uncertainty on the background yield is taken into account by varying the normalization and the 
shape of the distributions of several contributing processes. For both IT-i-jets and Z-i-jets production, the 
uncertainty on the normalization is studied following the recommendations in ref. [70]. For the 5-jet final 
state, a total uncertainty of 54% is assessed [70]. For the VT-i-jets background, the shape uncertainties 
due to the events with jets originating from heavy-flavour quarks are studied by varying the fraction of 
these events in the sample. The evaluation for this analysis follows the method described in ref. [40]. 
The shape and normalization uncertainties on the misidentified and non-prompt lepton component are 
propagated to the top-quark mass. The most important background topologies originate from single-top 
plus jets production. The impact on the top-quark mass is estimated by comparing the nominal yield 
(obtained using the Powheg generator interfaced to Pythia) with the equivalent result with a different set 
of generators (MC@NLO simulation for the s- and ITt-channels and AcerMC [71] for the generation of 
t-channel events). The effect of the (MC) top-quark mass used in the single-top background evaluation 
is also estimated by using two different input masses: 172.5 and 175 GeV, as well as differences in the 


16 


kinematics of the single top events. The uncertainty quoted as background in table 3 is calculated as the 
quadratic sum of all the above contributions. 

Table 3 summarizes all uncertainties on the estimated top-quark pole mass. 


Description 

Value 

[GeV] 

% 

pole 

nr 

t 

Statistical uncertainty 

173.71 

1.50 

0.9 

Scale variations 

(+0.93, -0.44) 

(+0.5, -0.3) 

Proton PDF (theory) and 

0.21 

0.1 

Total theory systematic uncertainty 

(+0.95, -0.49) 

(+0.5, -0.3) 

Jet energy scale (including Z^-jet energy scale) 

0.94 

0.5 

Jet energy resolution 

0.02 

<0.1 

Jet reconstruction efficiency 

0.05 

<0.1 

Z?-tagging efficiency and mistag rate 

0.17 

0.1 

Lepton uncertainties 

0.07 

<0.1 

Missing transverse momentum 

0.02 

0.1 

MC statistics 

0.13 

<0.1 

Signal MC generator 

0.28 

0.2 

Hadronization 

0.33 

0.2 

ISR/FSR 

0.72 

0.4 

Colour reconnection 

0.14 

<0.1 

Underlying event 

0.25 

0.1 

Proton PDF (experimental) 

0.54 

0.3 

Background 

0.20 

0.1 

Total experimental systematic uncertainty 

1.44 

0.8 

Total uncertainty 

(+2.29, -2.14) 

(+1.3, -1.2) 


Table 3; Value of the inferred top-quark pole mass and breakdown of their associated uncertainties. 


8.5 Studies on the definition of the extra jet 

The extra jet of the tf-t 1-jet system is required to have a pj larger than 50 GeV but other possibilities were 
also investigated. The full analysis was repeated with the pj of the extra jet satisfying different conditions 
such as pt > 30 GeV and pj > 40 GeV. The results differ from the baseline central value by less than 0.1 
GeV but the change in the systematic and statistical precision of the measurements was significant. As 
the pt requirement was decreased some systematics uncertainties, such as that due to the JES, increased 
and the statistical uncertainty became smaller. The original pT condition of 50 GeV represents a good 
compromise for the overall balance of these uncertainties and therefore was used as the baseline of the 
analysis. 
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8.6 Uncertainties on the measured “/^-distribution 


The experimental uncertainties of the unfolded “/^-distribution for each interval of ps are listed in table 4. 
These uncertainties are obtained following the same methodology as described in previous sections for 
the top-quark mass but this time applied to the ^/^-distribution. For each bin of the ^/^-distribution, all 
uncertainties have been added in quadrature. 


Ps interval 

n 

Stat. Unc. (%) 

Syst. Unc. (%) 

0 to 0.25 

0.126 

12.8 

7.1 

0.25 to 0.325 

1.122 

6.6 

4.5 

0.325 to 0.425 

2.049 

5.0 

3.5 

0.425 to 0.525 

2.622 

4.6 

2.1 

0.525 to 0.675 

2.125 

4.1 

3.1 

0.675 to 1.0 

0.302 

8.2 

8.1 


Table 4; Measured values of the !/?-distribution and their experimental uncertainties in percent. The statistical 
uncertainties are derived from the covariance matrix of eq. 9. 


9 Results and discussion 

This paper describes an experimental measurement of the top-quark mass using the novel method pro¬ 
posed in ref. [14]. The value of is obtained from a fit to the normalized differential cross section 
for tt production with at least one extra jet, tt + 1-jet, as a function of the inverse of the in¬ 
variant mass of the tt + 1-jet system, ps. This method allows a rigorous theoretical interpretation of the 
extracted mass parameter in terms of the top-quark pole mass or the running mass in the MS scheme. In 
the present analysis only the top-quark pole mass is measured, although future studies should also 

be able to determine the running mass when the theoretical calculations become available. 

Events with the tt + 1-jet final stale were selected using 4.6 fb“^ of 7 TeV pp collision data collected 
by the ATLAS experiment at the LHC in 2011. The total background in the tt + 1-jet sample is ~6%. 
Many distributions were studied to demonstrate the overall agreement between the MC predictions and 
data. A thorough study of the systematic effects with impact on the measurement was carried out and the 
associated uncertainties quantified. Experimental systematic uncertainties were computed for detector, 
signal and background modelling. 

The measured top-quark pole mass is: 


= 173.7 ± 1.5 (stat.) ± 1.4 (syst.)+^ ° (theory) GeV, 

where the theoretical uncertainties include the uncertainty due to missing higher orders in the perturb¬ 
ative NLO calculation, as well as uncertainties due to the PDL and as used in the calculations. The 
experimental uncertainty accounts for the uncertainties due to the imperfections in the modelling of the 
detector response, the background yield and the uncertainties arising from the signal modelling including 
hadronization. The dominant experimental uncertainties are due to the jet energy calibration (0.94 GeV) 
and the initial- and final- slate radiation modelling (0.72 GeV). 
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This measurement constitutes the first extraction of the top-quark pole mass from a measurement of the 
differential tt + 1-jet production cross section as a function of the inverse of the invariant mass of the 
tt + 1-jet system. It represents the most precise measurement of the top-quark pole mass to date with a 
total uncertainty of -1-2.3 GeV and -2.1 GeV. The value obtained for the top-quark pole mass agrees with 
the most accurate previous top-quark mass measurement in the pole-mass scheme [11] and with the direct 
top-quark mass measurement [6]. 

The analysis presented in this paper is statistically uncorrelated from the measurement using the 
inclusive cross-section measured in dilepton events [11]. The measurements could therefore potentially 
be combined, however this would require a detailed study of the correlations of both the uncertainties on 
the experimental measurements and the theoretical calculations. 
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